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Abstract. We introduce a measure of "quantumness" for any quantum state in a 
finite dimensional Hilbert space, based on the distance between the state and the 
convex set of classical states. The latter are defined as states that can be written as 
a convex sum of projectors onto coherent states. We derive general properties of this 
measure of non-classicality, and use it to identify for a given dimension of Hilbert space 
what are the "Queen of Quantum" states, i.e. the most non-classical quantum states. In 
three dimensions we obtain the Queen of Quantum state analytically and show that it 
is unique up to rotations. In up to 11-dimensional Hilbert spaces, we find the Queen of 
Quantum states numerically, and show that in terms of their Majorana representation 
they are highly symmetric bodies, which for dimensions 5 and 7 correspond to Platonic 
bodies. 
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1. Introduction 

The advent of quantum information theory has led to substantial efforts to understand 
the resources which are responsible for the enhanced information processing capabilities 
of quantum systems compared to classical ones. A large part of that research 
has been directed towards the creation and classification of entanglement [1, 2]. 
Entanglement plays an important role in quantum teleportation [3] and various quantum 
communication schemes [4] . It is also known that any pure state quantum computation 
which does not produce large scale entanglement can be simulated efficiently classically 
[5]. Physically, entanglement manifests itself as increased correlations between different 
subsystems compared to what is possible classically [6, 7]. But even for a system 
consisting of only a single subsystem one may ask how "quantum" a given state is, 
and what benefits one might draw from its "quantumness". 

In physics there is a wide consensus that the "least quantum" (or "most classical" ) 
pure states are coherent states. These are states which present the smallest possible 
amount of quantum fluctuations, as defined by a suitable Heisenberg uncertainty 
relation, evenly distributed over a pair of non-commuting variables. For example, in 
quantum optics, coherent states have that property of minimal and equal uncertainty 
for the field quadratures. Moreover, the dynamics of the latter is identical to that given 
by the classical equations of motion of the harmonic oscillator, and the property of 
minimal uncertainty is conserved during the time evolution created by the Hamiltonian 
of the electromagnetic field. The most classical mixed states possible can be obtained 
as a statistical mixture of coherent states. Any mixed state can be expanded over 
projectors on coherent states, with real coefficients given by the so-called P-function. 
Consequently, in quantum optics, states with a positive P-function are considered as 
"most classical" [8, 9]. 

In [10] we extended that definition to systems with Hilbert spaces of finite dimension 
d. These systems are formally equivalent to a spin of size j with d — 2j + 1. The P- 
function is not uniquely determined in that case, but leaves a lot of freedom in the 
specification of the higher spherical harmonics components. We therefore defined the 
set C of classical states (also called the set of "P-rep" states in [10], for positive P- 
function) as the ensemble of all density matrices p c for which a decomposition in terms 
of angular momentum coherent states \9ip) (see (8) for a precise definition) with positive 
weights exists, 



with < fa < 1 and X^/ 1 * = 1- At most d 2 = (2j + l) 2 terms in the convex sum are 
needed [10]. We showed that for j = 1/2 all states are classical in that sense. For spin 
j = 1 we found a necessary and sufficient criterion for classicality, and for higher values 
of j we found "non-classicality witnesses" which allow to easily detect large classes of 
non-classical states through the violation of necessary conditions for classicality derived 
from (1). For composite systems, the set of classical states C is in general strictly smaller 




(1) 
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than the set of separable states [10]. 

The definition (1) of a classical state allows so far just to determine whether a state 
is classical or not. However, it would be interesting to know how "non-classical" (or, 
in other words, how "quantum") a given state is, as one might expect that very non- 
classical states might be more useful for applications in quantum information processing 
than states which are only slightly non-classical. The situation is very analogous to the 
one encountered in the study of entanglement, where one wants to have an entanglement 
measure in addition to entanglement criteria. 

This is the question we are going to pursue in this work. We introduce a measure of 
quantumness in the next section, study some of its properties, and then apply it to find 
the "most quantum" states possible for a given Hilbert space dimension. We show that 
the largest possible quantumness can always be found in a pure state. The states with 
maximal quantumness turn out to possess remarkable geometrical beauty. We term 
them "Queens of Quantum" (QQ) states. In the lowest- dimensional non-trivial case 
(j — 1, i.e. d — 3) there is a unique Queen of Quantum (up to rotations of the coordinate 
system), which we determine analytically. In higher dimensions (up to j = 5) we find 
the Queens of Quantum numerically using quadratic optimization. Other attempts to 
define the "least classical" quantum states were proposed in the literature based on 
properties of the average value and the variance of the (pseudo-) angular momentum 
operator J [11, 12, 13]. We will briefly discuss these results in relation to the QQ states. 

2. Measure of Quantumness 

2.1. Definition and properties 

We define the "quantumness" Q(p) of an arbitrary state p as the distance from p to 
the convex set of classical states C. We thus introduce the measure of quantumness by 
defining 

Q(p) = mm ||p- p c | | , (2) 

PcSC 

where the minimum is over all classical states defined in (1), and = ti(A^ A) 1 / 2 
denotes the Hilbert-Schmidt norm. Note that our definition of quantumness is very 
analogous to the entanglement measure based on the distance of a state p from the 
convex set of separable states [14]. 

Several consequences follow immediately from (2): 

1. For any state p, and any dimension d, we have the bounds < Q(p) < 
i/trp 2 + A/trpf < 2. The lower bound is trivially realized for classical states. This 
implies Q(p) = for all states of a spin 1/2. 

2. An improved upper bound on Q(p) that only depends on the purity of p can 
be found by considering the distance to the maximally mixed state l/(2j + 1), which is 
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always classical [10]. This immediately leads to 

Q(p) < (3) 
For a pure state this bound coincides with the less stringent 

3. A different upper bound can be found by minimizing over a single pure coherent 
state \a) = \9ip): 

Q(p) < min||p- \a)(a\ \ \ < (1 + trp 2 - 2 max HJa)) 1/2 , (5) 

a a 

where H p (a) = (a\p\a) denotes the Husimi function of the state. For a pure state the 
bound becomes Q(\ip){ip\) < y/2(l - max a \{a\i))\ 2 ) 1 ' 2 . 

4. As the distance to a convex set is a convex function (see e.g. example 3.16 
of [15]), Q(p) is a convex function, i.e. for any two states p±, p 2 , and < p < 1, we have 

QiPh. + (i - p)P2) < pQ(pi) + (i - p)Q(p2) ■ (6) 

This implies that the quantumness of any mixed state cannot be larger than the largest 
quantumness of the pure states of which it is a mixture. 

5. Q(p) is invariant under rotations of the coordinate system. Indeed, let 
R n = exp(m.J) be a unitary operator associated with a rotation of the coordinate 
system about the axis n by an angle |n|. Since R n is unitary, we have \\p — p c \\ = 
1 1 RnpRn — RnPcRn \ \ for all density matrices p,p c - Furthermore, for p c e C, we also have 
p c = R^pcR^ G C, as R n only rotates coherent states into other coherent states, and 
therefore does not change the classicality of p c . Moreover, the map p c — > p c for given 
R n is an isomorphism C — > C. Therefore, we have min PceC = minp c£C . This leads to 
Q(p) = min Pc6C \\p - p c \ \ = min^c | \R n pR n - p c \ \ = Q(p) for p = R n pR n . 

With the same argument one shows that for a composite system consisting of s 
subsystems, Q(p) is invariant under independent rotations for all the subsystems, i.e. 
under transformations R = R ni £g> R n2 <8> • • • R ns ■ In addition, for a system consisting of 
s qubits, Q(p) is invariant under all local unitary transformations, as in this case local 
unitary SU(2) transformations leave the set C invariant. Our measure of classicality 
shares this property with any measure of entanglement. Indeed, for a multi spin- 
1/2 system, the set of classical states C is identical to the set of totally separable 
states. Therefore, in this case Q(p) coincides with an entanglement measure. For 
higher-dimensional subsystems, this is, of course, not true, as for example a SU(3) 
transformation of a single spin 1 can transform a coherent state into a non-classical 
state. 
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Figure 1. Quantumness of the state (7) as function of x = \plj sin#. The plot of (28) 
coincides with this curve within numerical accuracy. 



2.2. Simple examples 

2.2.1. Spin- 1/2 case As mentioned in [10], any spin-| pure state is classical, thus the 
set C coincides with the set of all quantum states. All states are thus trivially at distance 
from C. 

2.2.2. Pure spin-1 case In order to illustrate the behaviour of the quantumness in the 
simplest non-trivial case, let us consider the family of pure spin-1 states given by 

l^> = -rrr^(l 1 '- 1 ) + a; l 1 ' > + l 1 ' 1 )) ' (7) 

V x + 2 

where \jm), —j<rn<j, are eigenvectors of the spin angular momentum operator J z 
with eigenvalue m. Let x = \/2/ sin#. In figure 1 we plot the quantumness Q{\^x){i^x\) 
obtained numerically using the method described in section 3.4. We see that the largest 
quantumness is obtained at 6 = (or x = oo), corresponding to the state |1, 0). We will 
prove in section 3.2 that this state is indeed the spin-1 Queen of Quantum state and 
has Q = >/3/8. 

2.2.3. Thermal spin states Consider a density matrix given by the thermal state 
p = exp(—f3H)/ti(exp(—f3H)) with inverse temperature (3 = 1/ksT and Hamiltonian 
H . For — > all energy eigenstates are equally populated, so that p is the (properly 
normalized) identity matrix po corresponding to the maximally mixed state. Since 
Q(Po) = 0) we have the intuitively appealing result that for sufficiently large temperature 
all thermal states become classical. Indeed, there is a finite temperature where this 
happens, as there is a finite neighbourhood around po (with a finite radius in all 
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directions) where all states are classical [10]. For (3 — > oo the quantumness depends 
on the quantumness of the ground state (s) of H . 

We will illustrate what can happen with two examples for spin 1. In this case we 
have at our disposal a necessary and sufficient condition of classicality, which has been 
obtained in [10] (see Sect. 3.2 for a presentation of this criterion). 

1. For H = J z the ground state is the pure coherent state |1, —1) and thus it is 
classical. Our classicality criterion shows that in this case all thermal states are actually 
classical. 

2. For H = .11 , the ground state is the non-classical state |1, 0), and the classicality 
criterion shows that there is a critical temperature T c = 1//3 C = l/ln2 above which the 
quantumness disappears. 

Therefore, while the quantumness at low temperature depends on the model and 
in particular on the quantumness of low-lying states, classicality is a universal property 
for systems at thermal equilibrium for T — > oo. 

3. Queens of Quantum 

We now address the question of finding out which states are the "most quantum" states. 
This is a highly non-trivial question, as it requires to find a state that maximizes a 
quantity defined as a minimum over a convex set. We first introduce some definitions. 

3.1. Definitions and general properties 

We define the "Queen of Quantum" (QQ) states as those states pqq G Af for which 
Pqq = maxp 6 ^ Q(p), where Af is the set of all physical density operators acting on a 
Hilbert space of given finite dimension. In other words, the Queens of Quantum are the 
"most quantum" states for a given Hilbert space dimension. 

3.1.1. Pure versus mixed In order to identify the largest quantumness possible for a 
given Hilbert space dimension, we can restrict ourselves to pure states, as according 
to (6) the quantumness of a state p cannot be larger than the quantumness of the 
pure states from which it is a mixture. This does not immediately imply, however, 
that all QQ states are pure. Indeed, suppose that two pure states j^i) and \if 2 ) 
have the same quantumness, Q{\^i) = Q(\4>2) {ip2\) ■ Then, according to (6), 
QG#i)(^i| + (! - P)\^2){^2\) < for < p < 1. Equality is possible in 

principle, unless Q(p) is strictly convex, which, however, need not be the case. For 
example one might imagine that Af has a flat surface containing an d I^XV^I, 

in parallel to a flat surface of C containing the two closest mixed states for \ipi)(ipi\ and 
IV^) (^2 1- In this case, all states mixed from the two pure QQ states will have the same 
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maximum quantumness. The problem of degenerate quantumness of the QQ states is 
in fact generic, as we have seen that all states obtained by rotation of the coordinate 
system have the same quantumness. Nevertheless, we can start off by determining all 
pure QQ states, and then try to determine whether states can be mixed from these 
which would give the same maximum quantumness. For j — 1 we will show explicitly 
that all QQ states are pure. 

3.1.2. Majorana representation We now restrict ourselves to the case where p = \ip)(ip\ 
is a pure state. A spin-j coherent state can be expressed as [16] 

|0» = (j + m ) 2 ( sin 0' +m ( cos I)' m e-^\jm). (8) 

Note that with this definition the state \j — j) corresponds to 9 = 0. We take this 

state as the South pole, and \j j) as the North pole. The state \9ip) can be seen as the 

state obtained by rotating \j — j) in a direction specified on the sphere S 2 by an angle 

9 about the y-axis followed by an angle about the z-axis. This means that the mean 

value (9ip\J\9ip) of the angular momentum vector points in the direction given by the 

usual (i.e. counted from the North pole) polar angle 9 = n — 9 and azimuth if. We draw 

the attention of the reader to the fact that in the whole paper 9 will be counted from 

the South pole. This perhaps unusual convention is adopted here in order to simplify 

subsequent expressions for the Majorana polynomials. 

Any pure state \ip) can be represented by its overlap with coherent states 

7 , 1 
1 V- ( 2 3 



(WHj^E^J^c- (9) 

where ( is a complex number defined by 

C = e*nan^. (10) 

The scalar product (9), up to a prefactor independent of is a polynomial of degree 
at most 2j in the complex variable (. If ipj ^ 0, it is a polynomial of degree 2j, and has 
2j complex roots Ci> • • ■ (2j- This so-called Majorana polynomial [17, 18] reads 

M (o = E ( 2 / ) V;C* = ^ ft^ - c*) • en) 

fc=0 ^ ' k=l 

The inverse stereographic projection of these 2j roots define 2j points on the unit sphere 
S 2 through (10). If the degree of M(() is D < 2j, the Majorana representation is defined 
as consisting of D points associated with the D roots of the polynomial and 2j — D points 
at the North pole, and the prefactor tpj in (11) is replaced by ipD-j- This set of points is 
called the Majorana (or stellar) representation. It entirely characterizes the normalized 
state j^) up to a global phase. A nice feature of this representation is that all points 
rotate rigidly if undergoes a rotation R n = exp(in.J) of the reference frame. If the 
state |^) is a coherent state, \ip) = \9ifi), then the 2j roots of its Majorana polynomial 
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are all equal and given by (k = tan £z^ e * They correspond to 2j points which are 
antipodal to the vector (9ip\J\9ip). For instance, the state \j j) (corresponding to 9 — ir), 
which is a coherent state whose mean value is located at the North pole, is represented 
by 2j points at the South pole (once more we recall that we count 9 from the South pole). 

There is a natural interpretation of the Majorana representation in terms of ten- 
sor products of spin states. As a spin-j state, \ip) can be seen as a fully symmetrized 
direct product of 2j spins \. Each spin-| state is a coherent state of the form (8). As 
such it can be represented by a Majorana point on the Bloch sphere in the direction 
antipodal to the coherent state. The state \ip) corresponds to a symmetrization of 2j 
Bloch spheres and thus can be represented by a set of 2j points on a single sphere. The 
Majorana representation is useful in many contexts: calculation of the Berry phase for 
pure states [19], proof of Sylvester's theorem on Maxwell multipoles [20], investigation 
of states which maximize the variance (AJ) 2 = (J. J) — (J) (J) [11], or states such that 
this variance is uniform over the unit sphere [12] . 

Since the set C of classical states is invariant under rotation (as coherent states are 
just rotated to coherent states) the distance from \if>) to C is the same for any rigid 
rotation of the 2j points on the sphere. The problem of identifying pure QQ states for 
spin j reduces to identifying optimal distributions of 2j points on the unit sphere, up 
to global rotation. 

3.1.3. Pure state as eigenfunction of its classical neighbour Before proceeding to the 
identification of the most quantum states, we derive an important property of QQ states. 
As we pointed out, the largest Q(p) can always be reached with a pure state, p = 
The squared distance from a pure state to C is maximized by the state ipQQ such that 

Q 2 (\^QQ)(^QQ\) = maxmin (l - 2 (if; \p c \ if;) + trp 2 ) . (12) 

m pc&c 

The necessary condition on state ipQQ is stationarity of Q 2 (\if>)(if>\) with respect to 
variations of \if>) constrained by the condition = 1. The first variation of the 

functional I [if;] = 1 — 2 (ip\p c \ip) + 2E (ip\ip), equal to 

51 = -2(5if>\p c \if>) + 2E(5if>\if>) + c.c , (13) 

where 2E is the Lagrange multiplier, must be zero for all \5if>). It follows that 

Pc \iP) = E\iP). (14) 

Consequently, the wavefunction of a QQ state is an eigenfunction of (the density matrix 
of) its nearest classical state. We checked that all QQ states obtained in the following 
sections indeed do satisfy this property. 

3.2. Spin-1 case 

We now turn to the analytic investigation of the simplest nontrivial case of spin 1. 
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3.2.1. Pure states Let J a , a = x,y,z, be the 3x3 angular momentum matrices for 
j = 1. One can expand any density matrix p in the basis of J a and (J a J b + -hJa)/^, as 

p = -l 3 + -u.J + - £ (W ab --5 ab ) . (15) 

a,b=x,y,z 

Here u is a real vector, W a real symmetric matrix with trace 1, and 1 3 the 3x3 identity 
matrix. They are related to p by 

u a = tr(pJ a ), W ab = ti[p(J a J b + J b J a )} - 5 ab . (16) 

According to [10], p is classical if and only if the real symmetric 3x3 matrix Z with 
matrix elements 

Z a b = W ab - u a u b (17) 
is non-negative. 

Let be a fixed pure spin-1 state. Its Majorana representation consists of two 
points on the sphere. Since any state obtained by a global rotation of these two points is 
at the same distance from C, one can without loss of generality consider that these two 
points are specified by angles (0, 0) and (ir — 9,0). In the \jm) basis, the corresponding 
state is given by (7), with sin6> = y/2/x, since the Majorana polynomial of state \^ x ) 
is given by M(() = (1 + \/2x( + ( 2 )/Vx 2 + 2. When 9 varies from to ir/2, x takes 
values between \/2 and oo. Any spin-1 state can thus be brought to the canonical form 
(7) with x e [y/2, 

Lemma. The state defined by two antipodal points on the Majorana sphere is at 
distance a/3/8 from C. 

Proof. Without loss of generality we take the two points to be at the North and South 
poles. The corresponding state is |1, 0). Its parameters in the expansion (15) are u = 
and W =diag(l, 1, —1). A coherent state has parameters u = n and W ab = n a n b , where 
n is a three-dimensional unit vector which we parametrize as 

n = (sin 9 cos </?, sin 9 sin ip, cos 9) T . (18) 

Since C is the convex hull of the set of coherent states, any element p c (u, W) of C can 
be written as 

N 

u = ^A.rii, (19) 
i=i 

N 

W ab =^\(n i ) a (n i ) b , (20) 
i=i 

where (rij) , a = x,y,z, are the components of the vector rii and N is an integer. The 
distance | \p — p c \\ 2 with p = |1, 0) (1, 0| is 

-^AiAjni.nj + -^A.A^ni.rijl 2 - - X iJ2 w ob(ni)a(rn)b + j ■ (21) 

i,j i,j i a,b 
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(22) 



- W ab{ni)a{ni) b - ^2 W ab{nj)a{nj)b + 3). 

a,b a,b 

We parametrize vectors rii by angles 6>j and as in (18). After some trigonometric 
simplifications we obtain 

\\p ~ Pc|| 2 = g + - ^ XjXj ( - sin 2 Oi sin 2 9j cos2(y9j - ipj) 

2 sin 6i sin 9 j cos(ipi — <fj)(l + cos 9i cos 9j) (23) 

3 3 \ 

2 cos 9i cos 9j + -(2 — sin 9i — sin 9j) + - cos 2 9i cos 2 0j J . 

All terms in (24) are positive or can be expanded as a sum of terms of the form 
Xififii)) 2 '. Thus ||p — p c 1 1 2 > 3/8 for all states p c . This minimum is reached if 
all terms in (24) vanish, which implies the conditions 9 { = it/2 and 

Aj cos <fi = 0, Aj cos 2(/?j = (24) 



+ 
+ 



Aj sin </?j = 0, Aj sin 2<pj = 0. 



(25) 



These equations admit the solution Aj = 1/N and ipi = 2m /N for iV > 3, which 
corresponds to state 



(26) 



/ I 

4 



Pc 



\ 







Since p c G C, the minimum 3/8 is indeed reached. Thus the distance between |1, 0)(1, 0| 
and C is y/t/8. □ 

Theorem. The state |1,0) is the unique pure QQ state up to rotations. Its Majo- 
rana representation is given by a pair of antipodal points. 

Proof. Since |1,0) is at distance a/3/8 from C, it suffices to show that any other pure 
state is at distance smaller than a/3/8, for instance by explicitly exhibiting a classical 
state which is closer. We distinguish two cases. For x > y/E one can show, using the 
Z-criterion (17), that the state 

^ \ a(x) b(x) ^ 



Pc{X) 



4 

a(x) 



2 a \ x ) 
y b(x) a(x) \ 

with a(x) = xj[x 2 + 2) and b(x) = \j[x 2 + 2) is classical. The distance 



(27) 



/ 



x)(^x\ - Pc{x) 



3 / x 



8 V x 2 + 2 



(28) 
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is a strictly increasing function of x on [y/E, oo[, and thus any state \ip x ) with y/E < x < 
oo is at distance squared less than 3/8. 

For x < y/E, one can similarly show that the state given by (27) with a(x) = 
xj [x 2 + 2) and b(x) = A(x/ (x 2 + 2)) 2 — 1/4 is classical, and that its squared distance to 

111^X^1 - Pc {x)\? = ( " 2 ;ff+2)^ 12) ' (29) 

is a strictly increasing function of x on [a/2, y/E], thus bounded by its value at y/E, which 
is 3/32. Thus, the state with x = oo and the states obtained from it by rotations are the 
only pure states at distance a/3/8, all other being closer, which completes the proof. □ 

Numerical evidence shows that the distance given by (28) is precisely the distance 
between and C for all x (see figure 1). 

3.2.2. Mixed states The above Theorem shows that only states of the form R\1,0), 
where R is a rotation of the coordinate frame, are pure QQ states for spin 1. Since pure 
QQ states are at a distance Q 2 = 3/8 from C, any mixed QQ state has to be at least at 
the same distance. But, as explained in point 4. of section 2.1, convexity implies that 
mixed states can never be further away from C than the pure states they are composed 
of. Therefore, only mixtures p of pure QQ states that verify Q(p) 2 = 3/8 are candidates 
for mixed QQ states. 

Let p be a mixed QQ state. Its most general form of p is 

p = ^it>|l,0)<l,0|4 (30) 

i 

where Ri represents an arbitrary rotation of the coordinate frame, and the pi are positive 
and sum up to 1. The state 

l\ o o\ 

\ R\ (31) 

v° -J 

belongs to C and, by convexity of the norm, satisfies \\p — p c \\ < a/3/8. Since 
Q(p) = a/3/8, p c is indeed the classical state closest to p, which implies that p c 
has to be on the boundary of C. From (16), coordinates u and W are linear in p. 
Moreover, they transform respectively as a vector and a 2-tensor under rotations (see 
e.g. [21]). Using the fact that coordinates for state diag(l/4, 1/2, 1/4) read u = and 
W =diag(l/2, 1/2, 0), for p c we have u(p c ) = and 

/ k \ 

R\. (32) 



y^^PiRj 



\ 



i \ J 

The set C is the set of density matrices such that Z > 0. Since p c is on the boundary of 
C its matrix Z, which is equal to W(p c ), has a vanishing eigenvalue. Thus there exists 
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some vector n such that J2 ab W a bn a nb = 0. Using (32) one easily concludes that for all 
rotations, i?jn is equal to (0, 0, 1) T . Thus either all rotations are equal or they have the 
same rotation axis (0, 0, 1) T . In the latter case the state |1, 0) is invariant under Ri, and 
in both cases one concludes from (30) that p is a pure state. 
Thus 1 1,0) is the unique spin-1 QQ state up to rotation. 



3.3. Higher values of j 

Again we concentrate on pure states. The problem reduces to finding the maximum 
over all pure states \vp) of the minimum over p c e C of 

tr (| V) (VI - Pc) 2 = 1 - 2 J2 A *l M«i> I' + E X M(ai\a k )\ 2 , (33) 

i i,k 

where \a.j) = \0i(pi) are coherent states. The last term in (33) involves the overlap 
between coherent states 

\(a\a')\ 2 = co^^f^, (34) 

where j(a, a') is the angle between the two points corresponding to (a|J|o;) and 
(a' | J | a'). The other sum in (33) involves the overlap between |V) and the coherent 
states 

1 U=l L - U& 2 i=l 

where Q are the Majorana points corresponding to state |V) (for simplicity of notation, 
we assume that Vj 0) an d l( a Xi) is the angle between the point (a|J|a;) and the 
Majorana point Q. Since Q are the roots of the Majorana polynomial, whose coefficients 
depend on the components V* of |V), it is possible to show, using coefficient-root relations 
and normalization of |V), that 

with Ufc the kth symmetric polynomial of the Q (a — 1, o\ — ^ Q, a 2 = Yli<j dCjv)- 
For the lowest values of j one can express this quantity as a function of the angles 7 ifc 
between points Q and ( k . It is equal to 1 for j = 1/2. For j = 1, 3/2, we have 

— ^— w = 1 — 77T > sin — . (37) 

rr 2j cos 2 ^ 2? ^ 2 v ; 

For j = 2, 5/2, 



Jl 2j cos 2 ^ 2? ^ 2 



1 V"^ • 2 7ifc -2 Ti'fc' 

sin — — sin 



+ On(On -~\\ E 



2j(2j-1) ^ . 2 2 

all pairwise 
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with the last sum running over all ways of pairing 2j points into two distinct pairs. 
These formulae should easily generalize to general j. The whole expression (33) can 
thus be expressed as a function of terms of the form sin(7/2), which are equal to half 
the Euclidean distance between a pair of points separated by an angular distance 7. For 
instance in the case j = 1 the problem corresponds to finding 

■ 2 jjajM) ciT1 2 jjoiM) 



2 2 ^ , /c 



Our quest for pure QQ states thus amounts to finding an optimal arrangement of 
points on the sphere with two types of particles Q and ctj. The problem of arranging 
points on a sphere as evenly as possible has a long history. It was known by the ancient 
Greeks that 4,6,8,12 or 20 points could be arranged in a regular way. About 25 centuries 
later, by classifying all finite subgroups of the group of rotations in M 3 , it was proved 
that only five regular polyhedra exist. In the framework of electrostatics, one can define 
a generalized Coulomb potential between n point charges on the sphere as 

E^T (40) 

where dij is the distance between points i and j, and m a positive integer. The question 
of finding a configuration of points on the sphere that minimizes the potential (40) 
was first investigated by Thomson [22]. Similar questions appear in many fields, from 
crystallography to biology (see [23] and references therein). Our problem bears some 
similarity with such questions. However, our potential is more complicated (see e.g. (39)) 
and two kinds of "particles" are involved. Intuitively, for fixed Q, the minimization 
problem in (39) corresponds to finding coherent states ctj as remote as possible from the 
Q and from each other. 

Given the complexity of the minimax problem of the kind of (39) beyond the case 
j = 1, we choose a numerical approach. Many algorithms were devised to numerically 
obtain optimal configurations of points. Rather surprisingly it turns out that the optimal 
distribution does not necessarily coincide with regular polyhedra even in the case where 
these exist (see e.g. [24], where the distribution of point charges that minimizes Coulomb 
potential (40) is given up to 60 points). In the next subsections we will apply numerical 
techniques to identify QQ states for the smallest Hilbert dimensions. 



3.4- Numerical procedure 

The problem of finding the QQ states can be reformulated in terms of convex 
optimization, and even as an instance of quadratic programming. 



3.4-1- Quadratic programming For a fixed state p G A/", we represent the matrix p c G C 
minimizing the distance to C as a linear combination of coherent states whose directions 
densely and uniformly cover the unit sphere, p c = J2f =1 \ \(Xi) (eti\ with N large. Q 2 (p) 
then becomes a quadratic function of the coefficients Aj which has to be minimized 
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under the constraints \ > 0, ^2 =1 \ — 1- This is a problem of quadratic programming 
which can be solved by a variety of algorithms. Although the original linear combination 
contains several thousand coherent states, only a few of them enter the solution with 
coefficients significantly different from zero. In the search of QQ states, the result of the 
quadratic minimization was then numerically maximized by variation of the pure state 
p. The Majorana configurations were found numerically and then deformed to closest 
symmetrical figures under the condition that Q increased. It is probably superfluous to 
stress that maximization and minimization do not commute, so that the maximin and 
minimax of the squared distance \\p — p c \\ 2 are different. 

The optimization algorithm itself starts from N coherent states randomly 
distributed on a sphere, on which quadratic programming is performed. This yields 
an intermediate optimal state which is a combination of a relatively small number 
M ~ 5 — 25, M <C N, of coherent states. If some of these coherent states point in 
directions closer than a certain threshold, say, 2 degrees, they are replaced by a single 
coherent state with a cumulative weight directed along the weighted average direction. 
This step yields M' coherent states. Then N — M' new random coherent states are 
generated and the quadratic programming algorithm is run again starting from these 
iV coherent states (N — M' new and M' old). Iterating the process K times (typically 
K 1000 — 5000) with iV 100 yielded a 7-8 digit precision of the squared distance value. 

As expected from the remark of section 3.1.3, the numerically found wavefunction 
of a QQ state for given j coincides, within numerical accuracy, with an eigenvector of 
the density matrix of its nearest classical neighbours. In fact, the accuracy with which 
this property was fulfilled could serve as a measure of accuracy in the search of the 
maximin. An interesting point is that all QQ states that we obtained were invariably 
associated with the eigenvector of the density matrix of the classical state corresponding 
to its largest eigenvalue. 

3.5. Results 

We carried out a numerical search of the QQ states for j from 1/2 to 5. The resulting 
values of the distance are plotted in figure 2 together with the upper bound (4). The 
distance of the QQ state from C grows monotonically with the increase of j; an almost 
perfect fit is 

1-2/(2^ + 1). (41) 

The resulting arrangements of the Majorana points on the unit sphere that represent 
the QQ states, as well as the respective sets of coherent states constituting the nearest 
p c are given in Table 1. The numerically obtained values for the Majorana points and 
coherent states are listed in the Appendix (Tables Al and A2,A3). For the first few 
values of j it is possible to identify regular structures from these numerical results. They 
can be recognized as highly symmetric figures (see figures 3 and 4) in which 9 and tp 
are typically rational multiples of n. Table 1 gives the analytical expression that we 
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Figure 2. Squared distance of the Queen of Quantum states from the boundary of the 
classical domain C for 2j = 1 — 10 (dots), the empirical approximation 1 — 2/(2 j + 1) 
(dashed line) and the rigorous upper bound 1 — 1/(2 j + l)(full line). 



obtain for the QQ states in the | jm) basis if we identify the numerical results with these 
regular structures. We find in particular that for j = 2 and j = 3 (i.e. 4 and 6 Majorana 
points), the QQ states correspond to the Platonic bodies (tetrahedron and octahedron, 
respectively). For j = 4, where one would expect to find the next Platonic body, the 
cube, the symmetry of the Majorana configuration is lower (see figure 3). 

The number of local maxima of Q 2 with close values rapidly increases with the 
growth of j while the maxima themselves tend to become more and more shallow. This 
makes the search of the optimal configuration more and more difficult while the value 
of the maximin itself can still be reliably established. 

It is instructive to compare our results with the optimal distribution of identical 
point charges on a unit sphere which interact through the standard Coulomb repulsion, 
plotted e.g. in [22]. The analogy between the problems follows from the possibility 
to express the optimized quantity in terms of distances between the end points of the 
Majorana vectors and the vectors of the coherent states (see section 3.3). In the range 
j = 1 — 5 the symmetry of the Majorana configurations of the QQ states coincides with 
that of the equilibrium configurations of 2j charges on a sphere, the only exception 
being j = 4. Indeed, the optimal Coulomb configuration of 8 identical point charges is 
the square antiprism with symmetry D44. However that latter configuration gives only 
a local maximum of Q 2 equal to 0.76868 which is slightly less than the global maximum 
0.77108 realized in a configuration with lower symmetry D 2 d (Table 1). 

It is also interesting to compare our results to the anticoherent spin states 
introduced in [12]. These states are defined such that (J) = and (AJ n ) 2 = 
(n.Jn.J) — (n.J)(n.J) is uniform over the unit sphere, i.e. independent of n. Platonic 
states for j = 2,3,4,6 and 10 are shown to be anticoherent. In [13] it was shown that 
multiqubit states with diagonal spin covariance matrix and maximal variances of each 
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Figure 3. Queens of Quantum: The Majorana points for j — 1.3/2,2 (top), 
j = 5/2, 3, 7/2 (middle), j = 4, 9/2, 5 (bottom). See Table Al for a precise definition. 

spin component are optimal for reference frame alignment. This property is verified for 
our QQ states for j = 2, 3. QQ states with j = 2, 3 are therefore both "anticoherent" 
and optimal for reference frame alignment. 

4. Conclusion 

In summary, we have introduced the "quantumness" Q(p) for any finite-dimensional 
quantum state p. This quantity measures how quantum a state is. Q(p) is a real 
valued, positive, and convex function with a value between and 1 (more precisely, 
< Q(p) 2 < 1 — 1/d, where d is the dimension of the Hilbert space) that measures the 



Quantifying Quantumness and the Quest for Queens of Quantum 



17 



3 


Qqq 


Majorana points of \ip) on 
unit sphere 


QQ state 


Coherent state points of 
closest classical neigh- 

nniirs 


1 


3/8 


Two antipodal points 


7l(lM> + 


1,-1)) 


Equilateral triangle in 
equatorial plane 


3/2 


9/19 


Equilateral triangle in equa- 
torial plane 


75 (11' D + 


2' 2/ J 


Two points on poles, 
equilateral triangle in 
equatorial plane 


2 


16/27 


Tetrahedron 




'§|2, -2) 


Overturned tetrahedron 


5/2 


0.645914° 


Two points on poles, equi- 
lateral triangle in equatorial 
plane 


72 ( 1> t) + 


2' 2'/ 


Two parallel equilateral 
triangles symmetric on 
both sides of equatorial 
plane 


3 


347/486 


Octahedron 


-^(|3,2> + |3,-2» 


Cube 


7/2 


0.743138 6 


Two points on poles, reg- 
ular pentagon in equatorial 

Til T) & 


(2' 2) + 


Z _5\) 

2' 2// 


Two parallel regular pen- 
tagons 


4 


0.77108 


Four points in plane with 
line of symmetry; remain- 
ing four points obtained by 
improper 7r/2-rotation S4 

d hnnr q t r TTi m pfn t 1 1 1 1 /_> 
dUULllj oy lillilt-ljl y llllt; 


see Table Al 


Twelve coherent states 


9/2 


0.79676 


Three triangles in parallel 
planes, central rotated by n 


see Table Al 


Two points on poles, 
two triangles symmetric 
with respect to equato- 
rial plane and three dou- 

nlpfc on pnn^tnr 


5 


0.81664 


Two points on poles and 
two squares in parallel 
planes rotated with respect 
to each other by 7r/4 


see Table Al 


Four squares in parallel 
planes; the two South 
squares rotated by 7r/4 
with respect to those in 
North hemisphere 



Table 1. Queens of Quantum: The Majorana points of the QQ states, the QQ states 
in \jm) notation, and the set of coherent states of the nearest classical neighbour. 
a Minimum of (270286 + 61910 cos(2a;) + 58680 cos(4a;) + 855 cos(6a;) + 1530 cos(8a;) - 
45cos(10a;) - 51200 sin(x) + 25600 sin(3a;) - 5120sin(5x))/262144 
b Minimum of (68477212 + 10990343 cos(2a;) + 18268726 cos(4x) + 2030189 cos(6a;) + 
845124 cos(8a;) + 25319 cos(10a;) + 26474 cos( 12a;) - 91cos(14a;) - 4014080 sin(x) + 
2408448 sin(3x) - 802816 sin(5a;) + 114688 sin(7a;))/67108864 
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Figure 4. Classical states closest to the Queens of Quantum: the coherent state 
points {6, if) for j = 1,3/2,2 (top), j = 5/2,3,7/2 (middle), j = 4,9/2,5 (bottom). 
See Tables A2 and A3 for a precise definition. 

Hilbert-Schmidt distance of p to the convex set of classical states, defined as states with 
positive P-function [10]. We have shown that thermal states always become classical 
(Q(p) = 0) for temperatures larger than a critical temperature that depends on the 
dimension of the Hilbert space and the Hamiltonian, whereas the ground state of a 
system may or may not have non-zero quantumness. We used Q(p) in order to find 
the "Queens of Quantum" states, defined as the states with maximum quantumness for 
a given Hilbert space dimension. Maximum quantumness can always be reached for 
pure states, and we have demonstrated that the Queens of Quantum states correspond 
to beautiful, highly symmetric bodies when expressed in terms of their Majorana 
representation. For the two lowest dimensions that allow for Platonic bodies (j = 2 
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and 3, with 4 and 6 Majorana points, respectively), they are indeed the corresponding 
Platonic bodies (tetrahedron and octahedron, respectively). For j = 4, lowering the 
symmetry allows to obtain an even larger quantumness compared to the one for the 
corresponding Platonic body (the cube), and we have identified numerically all other 
Queens of Quantum states for 3/2 < j < 5 using quadratic programming. 
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Appendix: QQ states and their closest classical state for j < 5 

This appendix lists the numerical results obtained from the algorithms explained in 
section 3.4. The Majorana configurations were found numerically and then deformed to 
the closest symmetrical figures under the condition that Q was increased. Results are 
shown in Table Al. The coherent states shown in Tables A2,A3 were numerically found 
for these symmetrical configurations. Contrary to the Majorana points, the symmetry of 
the coherent states was not enforced. Along with the numerically found values, we also 
give the likely exact values for the coherent states and the weights in Tables A2,A3, as far 
as they can be deduced from the numerically found ones and symmetry considerations. 
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Table Al. Numerical coordinates (6*, ip) of the Majorana points for the Queens of 
Quantum for 1/2 < j < 5. The Queens of Quantum are given by (0, ip) through 
(10),(11). 
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Table A2. Numerical coordinates and weights (9i,ipi,\i) of the coherent states for 
the classical states achieving the minimum distance to the Queens of Quantum for 
1/2 < j < 7/2. Namely, p c = 'Ei \i\$m)($m\ with \0 m ) given by (8). The 
(Of,(pf,Xf) are the presumable exact values deduced from the numerical data, based 
on symmetry considerations. In various cases their entry reads n.a.e., which means 
that no analytical expression could be extracted from the numerical data. In the case 
j = 1, the exact analytical solution (26) is given. 
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